Ordering of dipolar Ising crystals 
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We study Ising systems of spins with dipolar interactions. We find a sim- 
ple approximate relation for the interaction energy between pairs of parallel 
lattice columns of spins running along the Ising spin direction. This relation 
provides insight into the relation between lattice geometry and the nature 
of the ordered state. It can be used to calculate ground state energies. We 
have also obtained ground state energies and ordering temperatures To from 
Monte Carlo simulations. Simple empirical relations, that give To for simple 
and body centered tetragonal lattices in terms of lattice parameters are also 
established. Finally, the nature of the ordered state and To are determined 

for Fes clusters, which crystallize on a triclinic lattice. 
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Dipolar interactions can lead to order at low temperature in magnetic |l| as well as in 
ferroelectric systems [@,[§. Interacting dipolar Ising systems (DIS) have recently become 
the subject of special interest because spin quantum tunneling has been observed in them 
It can take place at temperatures that are well below their ordering temperature 
To (see below). The type of order that ensues and To are not trivially determined owing 
to the long range nature of dipolar interactions and to their sign changes. Luttinger and 
Tisza were able to show long ago that the type of ordering depends on the geometry of 
the lattice |J. The theory was later generalized by Niemeijer and Blotte 0. It enables 
one to determine, through laborious calculations, the ground state energies and types of 



order for simple Bravais lattices with up to two identical magnetic dipoles per unit cell. 
Ordering temperatures are known from Monte Carlo (MC) simulations for simple cubic 
lattices || . Unfortunately, symmetry forbids the existence of cubic dipolar Ising systems in 
nature. There are however magnetic systems, such as the so called single-molecule magnets 
T2 \ as well as some rare earth sulfides [[J, which (1) crystallize in other lattice structures 



and have sizable single-axis anisotropies, and (2) the organic material around these clusters 
precludes all magnetic interactions except for the magnetic dipolar one. Consequently, they 
are expected to behave as DIS, and they indeed seem to do so. |pHTl|] Magnetic relaxation of 
some of these single-molecule magnet crystals (Mni2 and Fes) have recently been simulated 
in connection with spin quantum tunneling experiments [11], but their thermal equilibrium 
properties have not, as far as we know, been studied. 

The aim of this report is twofold. First, we wish to point out that there is a simple 
approximate relation for the interaction energy between pairs of parallel lattice columns of 
spins running along the Ising spin direction, and that the nature of the ordered state follows 
from this relation. Second, we report Monte Carlo results for the ordering temperatures 
of DIS on cubic lattices, simple and body centered tetragonal lattices (that include Mn 12 
cluster crystals as well as on the lattice (triclinic) on which Fe§ crystallizes ]T5[ . 



We treat here systems of spins, on lattices to be specified below, that are restricted 
to point along the z-axis, and interact among themselves only through magnetic dipole 
interactions. More specifically, let Si = ±S be a spin on the z-th site of a simple cubic or 
tetragonal lattice, with Hamiltonian 



y 

where s« = Si/S, Vij = (a 2 /r i: ,) 3 (l — 3z?-/r?-), a z is the nearest neighbor distance along the 
z-axis, Tij is the distance between sites % and j, 

v = ^sSf/al (2) 

g is the gyromagnetic ratio, /xg is the Bohr magneton, and fiQ = Att x 10~ 7 in SI units. 



^From here on, unless explicitely specified to be in Kelvin, energies and temperatures are 
given in terms of vq. It may be helpfull to keep in mind that, (/io/47r)/i^/(lA) 3 ~ 0.622K. 



We have used the standard Metropolis algorithm in our Monte Carlo simulations [II 
We have mostly used periodic boundary conditions (PBC). For PBC, a spin at site i is 
allowed to interact only with spins that lie within a system sized box centered on site i. We 
have also simulated some systems with free boundaries, some shaped like a box, and some 
bounded by spherical surfaces. With free boundaries, all spins in the system are allowed to 
interact. The results obtained are in agreement with the expected behavior that follows from 
Griffith's theorem: that the thermodynamic limit is independent of boundary conditions and 
of system shape if no external field is applied [ |15|j . 

We next find the field B\ that an infinitely long column of spins pointing up, as in Fig. 
1, produces a distance a x away from it. First note that, as can be easy checked, B\ would 
vanish if the moment density were uniformly distributed rather than on a lattice. An exact 
expression, as an expansion in powers of modified Bessel functions, can be obtained for B\ 
by writing the magnetic dipole linear density along a column as a Fourier series with period 
and phase that depend on a z and d z , respectively (see Fig.l). The zeroth order term (which 
corresponds to a constant magnetic dipole density) does not contribute. The leading term 
is of the form (a x / a z )~ l l 2 exp(— 2na x /a z ) cos(2nd z /a z ). We have checked numerically that 

0.1357 



B exp[27r(l - a x /a z )] cos(2ird z /a z ) } (3) 



where 



Bo - f (4) 
471 ai 



is within 1% of exactness for a x /a z > 0.8. (Deviations from exactness increase as a x /a z 
decreases, up to 11% for a x /a z = O.5.). A B x < value means that a column of spins up 
produces a field that points down. 

On the other hand, the field B s of an infinite column of spins, all pointing up, at one 
of its own sites is easily found to be B s = 4. 808. ..So- This is much larger than B\ for 



any reasonable distance between columns, and gives therefore rise to ferromagnetic order 
within each column of spins. As shown below, this picture holds for tetragonal lattices for 
a x /a z > 0.6, but may break down for smaller values of a x /a z . 

One therefore expects the ordered state on simple cubic and primitive tetragonal lattices, 
for which d z = in Eq. (^|), to be a two-dimensional anti- ferromagnetic array of ferromag- 
netically ordered columns spins as shown in Fig. lb. On the other hand body-centered cubic 
and tetragonal lattices, in which d z = a z /2 for nearest neighbor columns, are expected to 
order ferromagnetically. This is as originally predicted by Luttinger and Tisza ||. 

Making use of Eq. (|3]) and the value of B s , the ground state energy e = 
— {\/2){gjj, B S){B\ + B s ) can be easily calculated for a simple cubic or tetragonal lattice. 
It is, 

e ~ -2.404^0- 1.8 gti B SB x . (5) 

This result is shown in Fig. 2, with data points obtained from Monte Carlo simulations, as 
a function of the basal plane lattice constant a x . For ferromagnets the calculation of Eq is 
somewhat more involved and is not attempted here. This is because for lines of finite length 
L z , equation (^) is applicable only if \n(L z /a z ) 3> a x /a z . For \n(L z /a z ) < a x /a z , B\ < 0, 
independently of d z , which gives rise to magnetic domains in ferromagnets. Equation @ 
can be used to obtain the "wall energy" that counterbalances "magnetostatic" energies in 



domain size and ground state energy calculations Jl6| . 

Data points obtained from MC simulations for the ordering temperature of DIS on prim- 
itive tetragonal lattices are also shown in Fig. 0a. The equation 



T ~2.5v (a z /a x y- 7 (6) 

provides the fit shown in Fig. ^|a. 

For body centered tetragonal lattices, ferromagnetic order ensues, as expected from the 
fact that d z = a z /2 then. The ground state energy is not simply obtained then, since there 
are long range contributions. Data points exhibited in Fig. |||b, show that the ordering 



temperature follows a slightly different rule from Eq. (P). As shown in the table, the best 
fit to T is then T ~ 5.8 v (a z /a x ) 2 - . 

Thermal equilibrium results obtained by MC simulations of DIS on simple cubic lattices 
of (1) up to 16 x 16 x 16 spins with PBC and (2) 2109 and 4169 spins on spherical systems 
with free boundary conditions are shown in Figs. |^a, |^b, and ^|c. These results, together with 
results (not shown) for smaller systems, as well as for box shaped systems with free boundary 
conditions of different sizes, lead to ordering temperatures and ground state energies that 
are, within statistical errors, independent of shape and boundary conditions. Similar results 
have been obtained for simple and body centered tetragonal lattices, from which the data 
points for T shown in Figs. Qa and 0b, respectively, were obtained. 



Triclinic crystals of Fe§ clusters provide an interesting example [T3|]. Their lattice ge- 
ometry is exhibited in Fig. |j. Note that lattice sites on columns marked with a diamond, 
circle, square, etc. are displaced along the z-axis a distance d z ps a z /4 with respect to sites 
on columns marked with a circle, square, triangle, etc., respectively. Thus, while all columns 
on the same row on the right hand side of Fig. ^ interact antiferromagnetically, two columns 
on adjacent rows do not interact, and two rows on alternate rows (e.g., on rows marked with 
circles and triangles) interact ferromagnetically. This implies that any given row of columns 
in Fig. 4 orders antiferromagnetically, but there is near degeneracy (if only nearest nearest 
neighbor interactions are taken into account) as to how adjacent rows order with respect to 
each other. How rows order with respect to each other is quite likely determined by longer 
range interactions, but this question is beyond the scope of this report. To find out how the 
system actually orders, we turn to MC simulations. The ordered state is depicted in Fig. 
4 as follows: ♦, •, ■, A, and T stand for columns with spin up, while 0, °, D, A, and V 
stand for columns with spins down. This system orders at a temperature given in table I. 

It is worth remarking that Eq. ([|) is expected to be applicable whether dipoles are 
point like or extended, or whether they are magnetic or electric, except that the constant 
0.1357i?o is different in each case. It's validity only requires that the wave number of the 
leading fourier component of the dipolar density along the relevant columns be equal to 



the smallest reciprocal lattice vector for the columns. Furthermore, Eq. (|3|) is applicable 
to other crystal structures, such as hexagonal, not explicitely treated here. It follows, for 
instance, that DIS, whether magnetic or electric, order antiferromagnetically in a primitive 
hexagonal structure, but order ferromagnetically in a hexagonal closed packed structure. 

In summary, we have shown that Eqs. ([5]) and (|5[) provide simple approximate relations 
for the interaction energy between pairs of parallel lattice columns of spins running along 
the Ising spin direction, and that the nature of the ordered state follows from this relation. 
/Rrom this simple relation, we have obtained ground state energies. We have also obtained 
ground state energies and ordering temperatures from Monte Carlo simulations. A simple 
empirical relation, Eq. (|6]) gives results for the ordering temperatures of DIS on simple 
and body centered tetragonal lattices. These results, as well as results for Fe 8 clusters that 
crystallize on a triclinic lattice, are summarized on table I. (Substituting g ~ 2.0, S = 10, 
and lattice constants (see Fig. f|) of Feg cluster crystals, gives T ~ 170 mK. Similarly, 
T ~ 430 mK is obtained for Mn 12 acetate, which crystallizes in a primitive tetragonal 
lattice with a z = 12.391 and a x = 17.32A 0.) 
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TABLE I. Type of ordered state, ordering temperature To, and ground state energy eo, for 



dipolar Ising systems on 


various lattices. 


All energies 


and temperatures are 


given in terms of vq, 


where v Q = (fj, /^TT)(gfx B 


Sf/a% and ( M /4tt) M |/(iA) 3 


~ 0.622K. 




LATTICE 


ORDER 


T 


£o 


validity 


SC 1 ' 2 


AF 3 


2.50(5) 


-2.68(1) 




BCC 1 


F 


5.8(2) 


-4.0(1) 




FCC 1 


F 


11.3(3) 


-7.5(1) 




Primitive tetragonal 


AF 3 


2.5(a z /a x ) 1 - 7 


Eqs. (3-5) 


a x /a z > 0.6 


BC tetragonal 5 


F 


5.8(a z /a x ) 2 - 


Fig. 2b 


a x /a z > 0.6 


Fes (triclinic) 


AF 4 


1.9(1) 


-2.73(3) 





1. Type of order and Eq (but not T ) first given in Ref. Q 

2. To, from MC simulations, given in Ref. ||. 

3. AF order depicted in Fig. 1. 

4. AF order depicted in Fig. 4. 



5. Tq, from MC simulations, given in Ref. [17| for a x /a z = 0.5. 
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(a) (b) 




FIG. 1. (a) Parallel lines of spins, showing distances a z , a x , and d z . (b) Antiferromagnetic order 
in simple cubic lattice. 
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FIG. 2. (a) Semilog plot of the thermal equilibrium ordering temperature To and of the ground 
state energy eq versus a x /a z of the magnetic dipole model on a tetragonal lattice, where a x and a z 
are the nearest neighbor distances along the basal plane and c-axis respectively, o and • stand for 
temperature and energy, respectively. The continuous line for energy follows from Eq. (j5|). The 
straight line through the temperature data points is a fit; the equation for it is The results for 
a x /a z = 1-4 are obtained from columns of 32 x 32 spins on the basal plane by 128 spins along the 
c-axis; the result for a x /a y = 1.5 is from a column of 8 x 8 x 512 spins. All values of e and T are 
given in terms of vq, defined in Eq. |2[ (b) Same as in part "a" but for a dipole model on a BC 
tetragonal lattice. The straight line through the temperature data points is a fit; the equation for 
it is 5.8(a z /a x ) 2 - . 
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FIG. 3. (a) The thermal equilibrium staggered magnetization m s versus temperature, o and • 
stand for systems of 8 x 8 x 16 and 16 x 16 x 16 spins, respectively, on a SC lattice with PBC. 
o and ■ stands for spherically shaped systems of 2109 and 4169 spins, respectively, with free end 
boundary conditions, (b) Same as in part "a" but for the staggered susceptibility Xs- (c) Same as 
in parts "a" and "b", but for the specific heat. 
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FIG. 4. Lattice constants of Fes cluster crystals are shown on the left hand side. The lattice 
parameters are taken from Ref. [13]. On the right hand side, the lattice as seen looking down the 
6-axis. Anisotropy constrains all spins to be either up or down along the 6-axis for T < IK in this 
system. The 6-axis is therefore our z-axis. 
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